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Abstract 

We determine the degree of fine tuning needed in a generalised version of the 
NMSSM that follows from an underlying Z4 or R symmetry. We find that it is 
significantly less than is found in the MSSM or NMSSM and extends the range of 
Higgs mass that have acceptable fine tuning up to Higgs masses of ~ 130 GeV. 
For universal boundary conditions analogous to the CMSSM the phenomenology 
is rather MSSM like with the singlet states typically rather heavy. For more 
general boundary conditions the singlet states can be light, leading to interesting 
signatures at the LHC and direct detection experiments. 



1 Introduction 



The recent results from the ATLAS and CMS collaborations provide an indication, albeit 
inconclusive, that the Higgs may lie in the range 124-126 GeV (l}]8). A mass in this range 
has significant implications for supersymmetric extensions of the Standard Model (SM) 
capable of ameliorating the (little) hierarchy problem. In particular in the Constrained 
Minimal Supersymmetric Standard Model (CMSSM) the fine tuning needed to achieve 
this mass is large, requiring a cancellation between uncorrelated parameters of order 1 
part in 300. In the more general context of the MSSM one still requires 1% fine tuning 
even for an extremely low messenger scale of 10 TeV |9[j] 

To accommodate a heavier Higgs while avoiding very large fine tuning requires new 
structure. For the case this results from new physics at a scale higher than the su- 
persymmetry breaking scale one can perform a model independent analysis through 



the inclusion of higher dimension operators 10 11 . This analysis showed that the fine 
tuning is extremely sensitive to an operator that is most simply generated through the 
integration out of a massive singlet chiral superfield that couples to the two Higgs su- 
permultiplets of the MSSM. The effect of this operator is to allow a Higgs mass as high 
as 130 GeV without increasing the low fine tuning found in the MSSM. This applies even 
if the mass of the singlet states is as high as 3 TeV. 

This structure is not the one found in the usual Next-to-Minimal-Supersymmetric- 
Standard-Model (NMSSM, see e.g. [12] for a review) which assumes there is an un- 
derlying Z 3 symmetry that forbids the singlet mass term and an explicit /i term. The 
reason this structure has been favoured is that it explains why the singlet should be light 
and, after supersymmetry breaking, the singlet superfield scalar component acquires a 
vacuum expectation value (vev) that generates a [i term of order the supersymmetry 



breaking scale. Fine-tuning studies of the NMSSM have been performed in e.g. 13 -17 



However recently it has been realised that if, instead of the Z 3 symmetry, the NMSSM 
has a discrete R symmetry, 7L^ or Z^, then, after supersymmetry breaking, both the 
singlet mass and the [i term are generated but both are constrained to be of order the 



supersymmetry breaking mass 18 19 . Moreover the R symmetry has the advantage 
that it forbids the dangerous dimension 5 proton decay operator^] and does not have 
the domain wall problem |20] that is associated with the 1i 3 symmetry of the normal 
NMSSM0 To distinguish the discrete R symmetric NMSSM from the usual NMSSM we 



1 Note that this definition of fine tuning differs from ours in the choice of the measure and the fact 
that the parameters are taken to be low-scale parameters. 

2 These operators, allowed in the MSSM and the usual NMSSM, must be suppressed by a mass of 
order 10 7 times the Planck scale! 



3 Avoiding unacceptable domain walls requires an R symmetry even for the usual NMSSM 21 22 
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will denote it by the GNMSSM (23 



In this paper we study the fine tuning in the GNMSSM in detail, without requiring 
the mass of the new singlet states be larger than the mass of the MSSM SUSY states. 



Also, in contrast to 23 , we take into account the complete superpartner mass spectrum 
at one-loop as well as all one-loop contributions and the dominant two-loop contributions 
in the Higgs sector. There is a broad minimum of the fine tuning for a Higgs mass 
between the LEP bound and 125 GeV. We discuss the phenomenology associated with 
these low-fine tuned points and the implications for dark matter abundance and direct 
dark matter searches. 



2 The GNMSSM 

2.1 The superpotential 

The most general extension of the MSSM by a gauge singlet chiral superfield consistent 
with the SM gauge symmetry has a superpotential of the form 

W = W Yukawa +^S 3 + (fi + \S)H u H d + ZS+^ s S 2 (1) 

= Wnmssm + nH u H d + £S + ^n s S 2 (2) 

where VWukawa is the MSSM superpotential generating the SM Yukawa couplings and 
Wnmssm is the normal NMSSM with a Z 3 symmetry. Here and in what follows capital 
letters refer to superfields while small letters refer to the corresponding scalar compo- 
nent. One of the dimensionful parameters can be eliminated by a shift in the vev v s . 
We use this freedom to set the linear term in S in the superpotential to zero, £ = 0. 



Such a superpotential can arise from an underlying 1^ or 1^ symmetry 18 19 . Before 
SUSY breaking the superpotential is of the NMSSM form. However after supersym- 
metry breaking in a hidden sector with gravity mediation soft superpotential terms are 
generated but with a scale of order the supersymmetry breaking scale in the visible 
sector characterised by the gravitino mass, m^i%. With these the renormalisable terms 
of the superpotential take the form [19] 

W z r ~ Wnmssm + m 2 3/2 S + m 3/2 S 2 + m 3/2 H u H d , (3) 

yV Z R ~ Wnmssm + m\/ 2 S (4) 

where the ~ denotes that the dimensional terms are specified up to 0(1) coefficients. 
Clearly the lif case is equivalent to the GNMSSM. After eliminating the linear term in 
S the 7L% case gives a constrained version of the GNMSSM with fi s /fi = 2k,/ X. 
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Note that the SUSY breaking also breaks the discrete R symmetry but leaves the 
subgroup Z^, corresponding to the usual matter parity, unbroken. As a result the 
lightest supersymmetric particle, the LSP, is stable and a candidate for dark matter. 

2.2 Supersymmetry breaking 

The general soft SUSY breaking terms associated with the Higgs and singlet sectors are 
U so ft = m 2 s \s\ 2 + m 2 hu \h u \ 2 + m 2 hd \h d \ 2 

+ (bn h u h d + \A x sh u h d + ^kA k s 3 + ^b s s 2 + £ s s + h.c^j . (5) 

Note that the shift in the vev v s that is used to eliminate the linear term in the super- 
potential does not imply that the corresponding soft term £ s is zero as well. 

These terms and the soft breaking terms associated with the squarks, sleptons and 
gauginos depend on the details of the supersymmetry breaking sector. Here we will 
first consider the simplest case, the CGNMSSM, chosen in analogy with the well-known 
CMSSM, with a universal scalar and gaugino mass and all other soft terms propor- 
tional to their corresponding superpotential couplings. In order to compare to the usual 
NMSSM case we will subsequently relax the universality of scalar masses to allow the 
Higgs masses to differ. We furthermore allow for independent A\, A K in the second case. 

The independent supersymmetry breaking parameters of the CGNMSSM are mo, 
TO1/2, A , B and £ s where A and B are the constants of proportionality associated 
with the trilinear and bilinear terms respectively. These parameters are defined at the 
unification scale, Mx, and must be evaluated at low scales using the renormalisation 
group running. 

Taking into account the supersymmetric parameters as well, the CGNMSSM within 
this simple supersymmetry breaking scheme is specified by the following set of parame- 
ters fi, fi s , A, k, mo, mi/2, A , B and £ s (in the more general case we have m^ d , m^, m 2 , 
A\, A K in addition). Trading Bq, £ s and \x for v, tan (3 and v s via the EWSB conditions, 
there are eight (thirteen) parameters defining these models. 
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3 Fine- Tuning and dark matter abundance: analysis 
methods 



3.1 The fine tuning measure 



As introduced in 24 , 25 , a quantitative estimate of the the fine tuning with respect to 



a set of independent parameters, p, is given by 

c^ln?^ t) 3v 2 

A = maxAbsAj, A p = — = \— . (6) 

amp v A op 

where v is the EW scakfl The quantity A" 1 gives a measure of the accuracy to which 
independent parameters must be tuned to get the correct electroweak breaking scale. 
The parameters, p, correspond to the nine (fourteen) independent parameters discussed 
above plus the top Yukawa coupling^] all defined at the unification scale and chosen to 
be of mass dimension 2 where appropriate, e.g. fi 2 . 

3.2 SUSY particle spectrum constraints from the LHC 

A large part of the low fine tuned parameter space is not viable because of the constraints 
from the LHC. If all squarks are roughly degenerate, as is the case in models with a 
universal squark mass at the high scale, the bounds on the squarks and gluinos are very 
stringent. The precise bounds depend on the details of the sparticle spectrum. Bounds 



for some cases can e.g. be found online at 27, 28 



The strongest bound assumes a light neutralino and that all squarks and gluinos are 
degenerate, giving m s - = trig > 1400 GeV, while if one assumes that the squarks are 
twice as heavy as the gluino, the bound on the gluino mass changes to nig ~ 700 GeV. 
If these assumptions are violated the bounds can be significantly weaker, e.g. a single 
stop below 250 GeV is still allowed! 

To allow for these uncertainties we take the more conservative superpartner bounds 
to be m > 1200 GeV for the first two generation squarks and the gluino and for the 



chargino we apply the bound > 94 GeV 29 . LEP also constrains the masses of 
sleptons to be above (9(100) GeV. However these limits are always fulfilled when we 
apply the cut on the squark masses because we assume a unification of all squark and 
slepton masses at the GUT scale. 



v — u„ + v d where v UtC t are the up and down sector Higgs vacuum expectation values. Here we 
work in conventions in which v ~ 246 GeV. 

5 We use the modified definition for fine tuning |26| for the top- Yukawa coupling, appropriate for 
measured parameters. 
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3.3 Implementation in SARAH and SPheno 



To perform an exhaustive study of the GNMSSM we have used the public Mathematica 
package SARAH 30-32 . SARAH contains a model file for a general singlet extended MSSM 



(SMSSM) which is easily adaptable to the GNMSSM. Using this model file SARAH ana- 
lytically calculates all mass matrices, vertices as well as the two-loop Renormalization 
Group Equations (RGEs) and one-loop corrections to self-energies and tadpoles. The 
calculation of the loop corrections is performed in DR scheme and 't Hooft gauge. The 
results can be used to get the entire one-loop corrected mass spectrum based on the ap- 
proach first used in [33] . For a detailed discussion of the use of this method to calculate 
the one-loop mass spectrum in extensions of the MSSM we refer to 
SARAH also provides an interface to produce modules for SPheno 



34 



36 



and 35 



37 ; all derived, 



analytical expressions are exported to Fortran code which can be compiled together 
with the public SPheno version. This offers the possibility to automatically get a fully- 
fledged spectrum calculator which outputs the SUSY one-loop mass spectrum based on 
a two-loop evaluation of all parameters as well as the widths and branching ratios of 
all sparticles and Higgs fields. The resulting SPheno version writes, for the parameter 
point under consideration, input files for HiggsBounds |38| which can be used to check 
all existing collider constraints in the Higgs sector. We have extended SPheno by in- 
corporating routines to calculate the fine tuning. As cross check we have implemented 
the equivalent functions in a MSSM version written by SARAH and compared the results 
with the fine tuning calculation of SoftSUSY [39]. In addition, we have linked the known 
dominant two- loop corrections in the Higgs sector involving the strong coupling, a s , and 



the third generation Yukawa couplings 40 43 included in the public version of SPheno 



For a comparison of the fine tuning within the NMSSM we have used the SPheno 

which includes also the known two-loop corrections of the 



version presented in 34 



NMSSM 44 



and implemented the routines for the calculation of the fine tuning. 

To this 



The dark matter relic density has been calculated with Mi cr Omegas 45f|47 
end we used the option of SARAH to write model files for CalcHep 



48 which can also 



be used with MicrOmegas. To interface SPheno and MicrOmegas we used the SLHA+ 
functionality of CalcHep 



49 



i.e. the information about the numerical values of the 
current parameter space point is passed to MicrOmegas via the output file written by 
SPheno in the SUSY LesHouches format 
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For our scans we have used SSP which is optimized for parameter scans using the en- 
vironment provided by model implementations in SPheno and MicrOmegas based on the 
SARAH output. SSP includes also routines to run Markov Chain Monte Carlos (MCMC) 
in order to find points in parameter space with a very high likelihood according to some 
defined constraints. 
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4 Exploring the GNMSSM 



In the following we will present the results of our scans over the parameters of the 
GNMSSM. We will start with the universal case, the CGNMSSM, and subsequently 
relax the requirement of universality. We are particularly interested in regions which 
allow for a rather large Higgs mass. The largest Higgs masses can be achieved when 
the additional tree-level contribution to the Higgs mass is large, corresponding to large 
A, (which implies smallish k [12]) and small tan/3. We randomly scan over all the free 
parameters within this region. 

4.1 The CGNMSSM 

In Figure [T] we show the fine tuning as a function of the lightest Higgs mass for both 
the CMSSM and the CGNMSSM, with and without imposing the constraints on the 
superpartner mass spectrum from the LHC and constraints from the relic density. All 
points have been subjected to internal consistency (e.g. no tachyons, correct EWSB 
etc.). In the case of the CNMSSM it is known that the requirement of universality is 
highly restrictive: both A and k have to be rather small and the Higgs mass is even 
smaller than in the MSSM case. We therefore consider the NMSSM case only for more 
general boundary conditions in the next section. Figure [T] shows that if no experimental 
constraints are imposed, the smallest fine tuning is achieved for small Higgs masses, as 
naively expected. However, when the parameter space is subjected to the experimental 
bounds, the lowest fine tuned part of parameter space is not accessible, which is the well 
known "little hierarchy" problem. 

As may be seen from Figure [T]the little hierarchy problem is significantly alleviated in 
the CGNMSSM for Higgs masses above 118 GeV and a large section of parameter space 
with only a mild tuning remains. Interestingly, larger Higgs masses up to 125 GeV do as 
well in terms of fine tuning. Given the vastness of the CGNMSSM parameter space, we 
did an iterative scanning procedure, randomly scanning a large part of the parameter 
space and then zooming into the interesting looking regions in several steps. Therefore 
the points shown in Figure [T] are not smoothly distributed over all of parameter space 
and one should not interpret the density of points as the probability of finding a viable 
point with a given fine tuning and Higgs mass. 

Of particular interest for the GNMSSM phenomenology is the SUSY conserving 
singlet mass parameter /i s , which sets the overall mass scale for the singlet and the 
singlino. For the CGNMSSM we show the lightest Higgs mass as a function of this 
parameter in Figure [2] It can be seen that for Higgs masses above 120 GeV, the GUT 
scale parameter fi s is rather large, implying that the phenomenology of the CGNMSSM 
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Figure 1: Fine-tuning vs. the lightest Higgs mass in the CGNMSSM case with universal 
boundary conditions for the CMSSM (orange), and the CGNMSSM (blue). The first plot 
corresponds to the unconstrained case, the second plot takes into account the LHC bounds on 
particle masses with a cut on squark and gluino masses of 1.2 TeV and the third plot assumes 
an additional upper bound on the neutralino relic density. 

is rather 'CMSSM like', as also observed in [23]. The difference to the CMSSM case is 
of course a larger Higgs mass and smaller fine tuning. We show three benchmark points 
of the CGNMSSM in Table [TJ After SUSY cuts are taken into account, for the Higgs 
mass region of interest the lowest fine tuned points we find have a fine tuning of about 
35. 

Dark matter 

If we require that the lightest supersymmetric particle (LSP) constitutes the dark matter 
of our Universe, additional constraints apply. In this work we will require a neutralino 
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Figure 2: The lightest Higgs mass vs. the SUSY conserving singlet mass parameter fj, s 
for a representative part of parameter space with fine tuning below 100. The light (dark) blue 
points are before (after) SUSY cuts are taken into account. 



LSP [j and that the relic density does not exceed the 5a WMAP-7 51 upper bound of 
Vth? < 0.1298. While an under-abundance could always be compensated by the relic 
density of a multitude of other particles, an overabundance would require a deviation 
from the standard thermal history of the Universe (or at least a sufficiently low reheating 
temperature, such that the dark matter candidate never reaches thermal equilibrium). 
It has been shown that it is not possible in the MSSM to get the preferred relic density 



and a fine tuning below 100 for Higgs masses above 120 GeV 52 . When we require 
the correct relic density for the CGNMSSM in addition to the LHC and LEP limits, 
the lowest fine tuning we find is about 50. For a large part of parameter space the 
lightest neutralino is rather bino-like and dark matter is overproduced, similarly to the 
CMSSM case. We find that all points that pass the relic density requirement correspond 
to the stau coannihilation region. A region similar to the focus point of the MSSM with 
a large Higgsino fraction of the LSP doesn't show up. The reason is that the region 
of interest in the GNMSSM corresponds to small values of tan /3 while the focus point 



prefers moderate or large values 53 . The direct detection cross section is typically 
below 5 • 10~ 47 cm 2 , well below the sensitivity of current dark matter searches. 



Another interesting scenario would be a gravitino LSP, which could make the region of parameter 
space corresponding to a stau NLSP viable. 
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4.2 Generalised boundary conditions 

In this section we will relax the universality condition on the GUT scale parameters 
and allow the Higgs and singlet soft masses , rrvf^ an d m% as well as the trilinear 
parameters A\ and A K to vary independently. In Figure [3] we show the fine tuning as 
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Figure 3: Fine-tuning vs. the lightest Higgs mass for the generalised boundary conditions for 
the MSSM (orange), the NMSSM (red) and the GNMSSM (blue). The first plot corresponds 
to the unconstrained case, the second plot takes into account the LHC bounds on particle 
masses and the third plot assumes an additional upper bound on the neutralino relic density. 

a function of the lightest Higgs mass, this time for the GNMSSM as well as the MSSM 
and NMSSM. For the GNMSSM we find no significant improvement of the fine tuning 
with respect to the universal case. Nevertheless, for large Higgs masses, the fine tuning 
is still significantly lower than in the MSSM or NMSSM case. One interesting differ- 
ence compared to the universal case is the freedom to have small (in the NMSSM case 
zero) values for the singlet mass parameter Therefore the singlet states can be light, 
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leading to potentially interesting phenomenology. One interesting possibility is that the 
lightest Higgs is mainly singlet and the second lightest Higgs corresponds to the mainly 
MSSM-like state. In Figure [4] we show the second lightest Higgs mass vs. the fine tuning 
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Figure 4: The second lightest Higgs mass vs. the fine tuning in the case the lightest Higgs is 
mostly singlet and evades all experimental bounds. The imposed cuts are 1.2 TeV on squark 
and gluino masses and an upper bound on the relic density according to WMAP. 

for the case the lightest Higgs is mostly singlet and evades all experimental bounds. We 
find that a light singlet can be an advantage for the dark matter abundance in that 
it seems much easier to reduce the relic abundance. In our scans much more points 
showed up with low fine-tuning, the correct relic density and with Higgs and squark 
masses in the preferred ranges for the case where the lightest scalar has a large singlet 
component. The origin of this is most likely the much weaker correlation between the 
singlino mass and annihilation cross section and the other masses and electroweak pa- 
rameters. Correspondingly we are no longer restricted to the narrow stau coannihilation 
region. Another interesting observation is that the direct detection cross-section is close 
to current bounds and hence this part of GNMSSM parameter space will be probed not 
only by the LHC but also by the next generation of direct detection experiments. We 
show two benchmark points for the generalised boundary conditions in Table [T] 

4.3 Comparison between the NMSSM and the GNMSSM 

It is interesting to explore the effect of the additional independent parameters of the 
GNMSSM compared to the NMSSM for the Higgs masses. To do this we have chosen 
to analyse a NMSSM point discussed in Ref. [54] with a Higgs mass near 126 GeV and 
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perturb away from the NMSSM point via the additional GNMSSM parameters. The 
results for a variation of (//, b/j.) and (/jl s ,/jl) are depicted in Figure [5j One may see that 
even small values of the new parameters can change the Higgs masses significantly Since 
the largest contributions of the fine tuning are due to the CMSSM parameters mo, mi/2 
and A , the variation of the parameters in the Higgs sector has only a very small impact 
on the overall fine tuning. Therefore, it is possible to vary the Higgs masses using these 
parameter without increasing the fine tuning. 



5 Phenomenological aspects of the GNMSSM 

The extension of the MSSM to include a singlet state can significantly change the phe- 
nomenology. The most important differences come about due to the possibility that 
there is an additional light Higgs state that is mainly singlet and the possibility that 
the LSP is mainly composed of the singlino. For the case of universal scalar boundary 
conditions the lowest fine tuned points correspond to large fi s (c.f. Figure [2j and in this 
case the singlet states are heavy and the phenomenology is close to MSSM phenomenol- 
ogy but with a Higgs that can be much heavier. However, for the case of generalised 



boundary conditions discussed in Section 4.2, low-fine tuned points can have small /i 



and hence different phenomenology while still having a heavy Higgs state consistent 
with the LHC hints for a 125 GeV Higgs. 



5.1 Higgs phenomenology 

Let us first consider the effect of a light, mainly singlet, Higgs, hi, on Higgs phenomenol- 
ogy. It can change the decay properties of the MSSM-like Higgs field dramatically as 
is depicted in Figure [6j Depending on the values of fi and fi a the decay of the mainly 
doublet Higgs, h 2 , into two singlet Higgs fields can be kinematically allowed. If this is 
the case the branching ratio Br(/i2 — > hihi) will be very large and the total width of the 
MSSM Higgs will typically be increased by more than one order of magnitude. Another 
interesting feature is that the up-type fraction of the MSSM-like Higgs is particularly 
sensitive to the [i parameter without significantly changing the down-type fraction. This 
allows for an increase Br(/i 2 — f 77) as may be indicated by the LHC measurements. 
As shown in the left plot of Figures [7j changing /1 from -20.9 GeV to -22.5 GeV (for 
/is = —55 GeV) enhances the branching ratio of the di-photon channel from less than 
2.6 • 10~ 3 to more than 2.9 • 10~ 3 . The reason is that the up-type fraction of the doublet 
Higgs changes by more than 10%. 

Higgs phenomenology can also be affected by the change allowed in the GNMSSM in 
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Figure 5: The mass of the lightest, singlet-like Higgs (left) and second lightest, MSSM-like 
Higgs (right) for a variation of the GNMSSM specific parameters (fi, b^) (first row) and (/i s , /i) 
(second row). The other parameters have been chosen as in [54]: mo = 385 GeV, m 1( / 2 = 
430 GeV, tan/3 = 3.0, A = 1590 GeV, A = 1.23, k = 1.07, A x = 1580 GeV, A K = 1560 GeV 
and v s = 250 GeV and the GNMSSM parameters not shown in the different Figures were 
kept 0. 

the neutralino structure, particularly in the nature of the LSP. As may be seen in Figure[8] 
there is minimal correlation between the singlino fraction or the singlino mass and the 
fine tuning. As a result it is possible to find areas in parameter space with low fine 
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Figure 6: Decay properties of the MSSM-like Higgs for a variation of [i and \x s . Left: Br(/i2 — > 
h\h\) (left), right: total width of /12 (right) in MeV. The other parameters were set to mo = 
640 GeV, m 1/2 = 480 GeV, tan/3 = 2.53, A = -253 GeV, A = 0.98, k = 1.19, A\ = 
136.2 GeV, A K = -410.3 GeV, v s = 620.8 GeV, by, = 95.2 GeV 2 , b s = 215.2 GeV 2 and 
£ s = —106.5 GeV 3 . The overall fine tuning for the shown parameter range is ~ 58 - 60. 
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Figure 7: Photonic decay of the MSSM-like Higgs. The left plot shows Br(/i 2 ->■ 77) x 1000. 
The other two figures show the up-type (middle) and down-type (right) fraction of the light 
MSSM-like Higgs. The parameters are the same as for Figure [6j Comparing to Figure [6] we 
see that the region where the photonic decay becomes very suppressed corresponds exactly to 
the region where the decay channel hi — > h\h\ opens up. 

tuning and m-LSP < \rfih- in this case the invisible decay h — > X1X1 is possible. Again, 
as shown in Figure [9| the exact branching ratio is also very sensitive to the parameters, 
particularly the value of fi s . A variation of /i s in the range [-50,50] GeV changes the 
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Figure 8: Singlino fraction of the LSP vs. fine tuning 




Figure 9: Decay of the light Higgs particle into two singlinos as a function of fi s . The left plot 
shows the mass of the LSP as function of /j, s , whereas the right plot gives the branching ratio 
Br(/i — > The other parameters have been chosen as mo = 525 GeV, m l / 2 = 591 GeV, 

tan/3 = 2.9, A = 1410 GeV, A = 1.42, k = [0.04., 0.05], A\ = -1815 GeV, A K = -569 GeV, 
v s = -360 GeV, \i = 60 GeV, bfi = 4056 GeV 2 , b s = -6545 GeV 2 and £ s = -144174 GeV 3 . 



mass of the MSSM-like Higgs only by 0.5 GeV, but can increase the branching ratio into 
two singlinos from less than 60% to nearly 100%. This is mainly caused by the increase 
in phase space due to the decrease in the mass of the singlino. If the indications of a 125 
GeV Higgs turned out to be true, these Higgs decays to lighter singlet states would be 
strongly constrained and the likely case that they are kinematically forbidden would be 
preferred. On the other hand, if the hints turn out to be a statistical fluctuation, these 
decays would offer interesting escape routes for Higgs physics. A more detailed study 
of these and other interesting aspects of the Higgs sector in the GNMSSM will be given 
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elsewhere 
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5.2 SUSY phenomenology 

The nature of supersymmetric signals is largely determined by the LSP. For the case 
that fx s is large the LSP is as in the MSSM so the signatures will be the same as in the 
MSSM. For the case that fi s is small and the LSP is mainly singlino the signatures can 
change. 

The LHC is sensitive to gluino and squark pair production. Let us consider first 
the case the gluino is lighter than the squarks. The dominant gluino decay mode is 
likely to be the three body decay into a quark anti-quark pair and the LSP. The singlino 
has no direct coupling to squarks so this decay will proceed via its neutral gaugino or 
Higgsino component. Over a large part of parameter space the singlino dominated LSP 
has a significant Higgsino component so the dominant decay mode will be into the third 
generation of quarks, g — > tis, bbs. Note that this expectation is different from the 
MSSM because the dark matter constraints on the makeup of the LSP are different. In 
the MSSM the LSP cannot be dominantly Higgsino and the gaugino component (mainly 
Bino) dominates. In the GNMSSM the singlino dominates but the Higgsino component 
is typically much larger than the gaugino component. Thus the collider signal for gluino 
pair production will be four jets associated with the third generation of quarks plus 
missing energy and momentum. The dominant decay of squarks is first to a quark and 
the gluino followed by the gluino decay into third generation quarks. Thus the signal for 
squark pair production will be multijets mostly associated with third generation quarks 
plus missing energy and momentum. 

For the case the squarks are lighter than the gluino their dominant decay mode 
will be the two body decay into the associated quark and singlino. Thus the signal for 
squark pair production in this case will be di-jet production plus missing energy and 
momentum. Gluinos will dominantly decay to a quark (anti) squark pair giving a final 
state of quark antiquark plus LSP. Thus the signal for gluino pair production in this 
case are fourquark-jets plus missing energy and momentum. 

6 Summary and Conclusions 

The realisation that simple discrete R symmetries naturally constrain the n and n s pa- 
rameters to be of the SUSY breaking scale mean that the GNMSSM is as natural a 
theory as the NMSSM. Indeed the fact that the symmetry also eliminates the dangerous 
dimension three, four and five baryon- and lepton-number violating terms in the La- 
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grangian and avoids destabilising tadpoles and domain wall problems renders it a more 
promising starting point than the NMSSM. 

In this paper we have shown that the GNMSSM, with gravity mediated SUSY break- 
ing and universal boundary conditions, can accommodate a Higgs mass up to 130 GeV 
without a significant increase in the lowest fine tuning needed which remains about 1 
part in 35. This is in contrast to the MSSM which requires a fine tuning of greater than 
1 part in 300 to accommodate such a Higgs mass. In this work we have assumed an 
underlying GUT structure at a scale around 10 16 GeV. In particular we have assumed 
universal squark masses mo at the high scale. If we were to relax this assumption, even 
smaller values of the fine tuning should be achievable, both due to a potentially lighter 
stop as well as due to a lower scale at which the fine tuning is evaluated, cf. e.g. (9]. We 
have also determined the dark matter abundance and shown that a significant propor- 
tion of the low fine tuned parameter space can lead to dark matter close to the WMAP 
bound without a significant increase in the fine tuning. The dark matter lies in the stau 
co-annihilation region and its direct detection cross section is far below the sensitivity 
of direct detection dark matter experiments. The low fine tuned points correspond to 
large fi s . In this region the new singlet states are heavy and the phenomenology is very 
close to that of the MSSM but with a heavier Higgs allowed. 

In order to make a comparison with the NMSSM we have also explored the more gen- 
eral case relaxing the universality condition on the GUT scale parameters and allowing 
the Higgs and singlet soft masses m^ d , m 2 hu and m 2 s as well as the trilinear parameters 
A\ and A K to vary independently. We found the fine tuning in the GNMSSM is rela- 
tively insensitive to this change and that it is significantly lower than that found for the 
NMSSM, particularly in the 125 GeV Higgs mass range. Interestingly for these more 
general boundary conditions the fine tuning is relatively insensitive to a reduction in the 
singlet fi s term and as a consequence there may be an additional light, mainly singlet, 
Higgs scalar and the LSP can have a large singlet component. This can change the 
phenomenology dramatically. The mainly singlet Higgs state can be the lightest state 
without conflicting with LEP and LHC bounds and indeed the doublet Higgs state can 
dominantly decay into pairs of singlet Higgs. Alternatively the doublet Higgs state can 
dominantly decay invisibly to a pair of mainly singlino LSPs. The case that the LSP is 
mainly singlino also affects the dominant decay channels of the gluino and squarks giv- 
ing characteristic signals that will provide tests of the scheme. Also a light singlino LSP 
changes the dark matter expectation as it is much easier to reduce the relic abundance 
so that we are no longer restricted to the narrow stau coannihilation region. In this case 
the direct detection cross-section can be close to current bounds and hence this part of 
GNMSSM parameter space will be probed not only by the LHC but also by the next 
generation of direct detection experiments. 
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BP1 


BP2 


BP3 


BP4 


BP5 


m [GeV] 


746 


163 


957 


573 


752 


mi/2 [GeV] 


476 


568 


557 


482 


472 


tan/3 


2.7 


2.9 


2.8 


3.4 


2.8 


A [GeV] 


1433 


1666 


782 


27 


-198 


A 


1.43 


1.47 


1.58 


1.34 


1.12 




-0.1 


0.09 


-0.005 


1.52 


1.03 


A A [GeV] 


A 


A 


A 


400 


192 


A K [GeV] 


A 


A 


A 


-323 


-326 


^ s [GeV] 


-841 


-190 


-929 


390 


281 


Us [GeV] 
™L [GeV 2 ] 
< [GeV 2 ] 
m 2 [GeV 2 ] 

/i [GeV] 
b\i [GeV 2 ] 

& s [GeV 2 ] 


-5931 


-5354 


-5799 


131 


-37 


m 2 
m 2 
m 2 
-750 


ml 
ml 
ml 
-1136 


ml 
ml 
ml 
-934 


9.1 • 10 5 
2.3 • 10 6 
2.8 • 10 6 
-33 


5.4 • 10 5 
2.4 • 10 6 
1.7- 10 6 
10 


-2.4- 10 6 


-1.2- 10 6 


-2.3 • 10 6 


147 


26 


— 1.9 • 10 7 


-5.4- 10 6 


-1.4- 10 7 


326 


144 


C [GeV 3 ] 


2.2 • 10 9 


1.5 • 10 9 


3.0 • 10 9 


22 


-8 


^squark [GeV] 


1256-1293 


1207-1263 


1507-1548 


1211-1248 


1280-1315 


m 5 [GeV] 


1219 


1389 


1416 


1242 


1235 


m hl [GeV] 


124 


123.5 


125 


93.5 


78 


m fe2 [GeV] 


1002 


856 


1257 


125 


124 


/ii singletfraction 


C(10~ 4 ) 


C(10~ 6 ) 


e>(io- 4 ) 


0.8 


0.85 


Br(h — >• 77) 


2.29 • IO- 3 


2.28 • 10~ 3 


2.2 • 10~ 3 


2.5 • lO" 3 


2.66 • lO" 3 


Br(6 ->■ S7) 


3.1 • 10~ 4 


3.1 • 10~ 4 


3.1 • 10~ 4 


3.1 • 10~ 4 


3.3 • 10~ 4 


Aa M 


-7.8- 10- 11 


-2.5 • 10- 10 


-5.4- 10- 11 


1.7- 10" 10 


8- lO" 11 


5p 


6.2 • 10~ 5 


6.6- 10~ 5 


7.5 • 10- 5 


1.9 • 10~ 4 


3.1 • HT 4 


m^o [GeV] 


229 


270 


168 


99 


70 


X? singlinofraction 


C(10~ 5 ) 


e>(io- 5 ) 


e>(io- 5 ) 


0.1 


0.2 


fi/1 2 


7.5 


0.10 


7.4 


0.017 


0.11 


a p [cm 2 } 


2.8 • 10" 47 


2.2 • 10~ 47 


6 • lO" 47 


1.2 • 10" 44 


1.3 • 10" 45 


A (Fine-tuning) 


34.9 


51.0 


51.8 


44.9 


52.7 



Table 1: Benchmark scenarios for the GNMSSM for the universal (BP1-BP3) and the general 
(BP4-BP5) case. m squar k shows the range of squark masses of the first two generations. For 
the last two points the second lightest Higgs is mostly MSSM-like. All input parameters except 
tan (3 and v s are given at the GUT scale. 
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